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ABSTRACT 

We report on the analysis of new and previously published MMT optical spectra of the black hole 
binary XTE J1118+480 during the decline from the 2000 outburst to true quiescence. From cross- 
correlation with template stars, we measure the radial velocity of the secondary to derive a new 
spectroscopic ephemeris. The observations acquired during approach to quiescence confirm the earlier 
reported modulation in the centroid of the double-peaked Ho emission line. Additionally, our data 
combined with the results presented by Zurita et al. (2002) provide support for a modulation with 
a periodicity in agreement with the expected precession period of the accretion disk of ~ 52 day. 
Doppler images during the decline phase of the Ho emission line show evidence for a hotspot and 
emission from the gas stream: the hotspot is observed to vary its position, which may be due to the 
precession of the disk. The data available during quiescence show that the centroid of the Ho emission 
line is offset by about —100 km s _1 from the systemic velocity which suggests that the disk continues 
to precess. A Ha tomogram reveals emission from near the donor star after subtraction of the ring- 
like contribution from the accretion disk which we attribute to chromospheric emission. No hotspot is 
present suggesting that accretion from the secondary has stopped (or decreased significantly) during 
quiescence. Finally, a comparison is made with the black hole XRN GRO J0422+32: we show that 
the Ho profile of this system also exhibits a behaviour consistent with a precessing disk. 
Subject headings: accretion, accretion disks — binaries: close — stars: individual: XTE J1118+480 
— X-rays: stars 



1. INTRODUCTION 

Low-Mass X-ray Binaries (LMXBs) are binary systems 
where mass from a low mass secondary (M2 < 1 Mq) 
is accreted onto a neutron star or black hole primary 
via an accretion disk. A subclass of LMXBs undergo 
episodic outbursts, separated by periods of relative low 
X-ray brightness, during which the X-ray luminosity can 
increase as much as a factor 10 7 . Known as X-ray Novae 
(XRNe) , these LMXBs have provided us with some of the 
most secure dynamically confirmed stellar black holes. 
For a review of the observational and theoretical status 
of LMXBs, and in particular XRNe, see van Paradijs 
& McClintock (1995), Tanaka & Shibazaki (1996), van 
Paradijs (1998) and McClintock & Remillard (2003). 

XTE J1118+480 is the first identified black hole XRN 
located in the Galactic halo (Wagner et al. 2001; Mirabel 
et al. 2001). After its discovery during a faint (39 
mCrab) X-ray outburst on March 2000 (Remillard et 
al. 2000), its optical counterpart was quickly identi- 
fied with a 12.9 magnitude emission-line object which 
had brightened by ~ 6 magnitudes from its pre-outburst 
value (Uemura et al. 2000a; Garcia et al. 2000). The 
multi-wavelength outburst spectrum of XTE J1118+480 
has been intensively studied (Hynes et al. 2000; McClin- 
tock et al. 2001a; Frontera et al. 2003). In particular, 
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thanks to the location of XTE J1118+480 close to the 
Lockman hole it was possible to carry out observations 
in the extreme UV (Hynes et al. 2000), a wavelength 
range inaccessible for the low galactic latitude XRNe. 

Optical photometry performed during outburst re- 
vealed a periodicity of 0.17078 d (Uemura et al. 2000b) 
which was later shown to be the superhump period 
caused by the precession of an eccentric accretion disk. 
During the decay of the outburst, McClintock et al. 
(2001b) and Wagner et al. (2001) measured the radial 
velocity curve of the companion star yielding to a mass 
function of ~ 6 M Q . This is beyond the maximum al- 
lowed mass for a neutron star (Rhodes & Ruffini 1974), 
a result which establishes the dynamical evidence for a 
black hole primary in XTE J1118+480. They also found 
that the XRN is at a distance of 1.8 ± 0.6 kpc and has a 
high inclination (i ~ 80°) relative to the observer. Sub- 
sequent spectroscopy has refined the classification of the 
secondary to a probable metal-poor K5-7 V subdwarf 
(McClintock et al. 2003). Additionally, the outburst UV 
emission line spectrum showed evidence for a nuclear- 
evolved donor star (Haswell et al. 2002). The mean 
rotational broadening of the secondary yields a mass ra- 
tio ofq = M 2 /M 1 = 0.037 ± 0.007 (Orosz 2001) which 
implies that the 3:2 resonant radius is within the disk 
tidal radius (see Whitehurst & King 1991 and references 
therein). Near quiescence, the optical light curves dis- 
played features characteristic of the ellipsoidal variations 
of the secondary star (McClintock et al. 2001b; Wagner 
et al. 2001). A 12 minute offset between the spectro- 
scopic and photometric ephemerides was also observed 
and tentatively explained by Wagner et al. as due to an 
asymmetrically emitting disk component. This hypoth- 
esis was confirmed by the observations of Zurita et al. 
(2002) who found a superhump component with a pe- 
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riod 0.3 per cent longer than the orbital period in their 
optical light curves obtained near quiescence. This su- 
perhump period implies a 52 day precession period for 
the accretion disk. 

In a previous paper, Torres et al. (2002a) presented the 
results of the analysis of their multi-epoch spectroscopy 
following the 2000 March outburst. They reported varia- 
tions in the double-peaked emission lines on a time-scale 
longer than the orbital period of the system. This was 
interpreted as evidence for the precession of an ellipti- 
cal accretion disk. During the approach of the system 
to quiescence, Zurita et al. (2002) found asymmetries in 
the Ha line and long-term variations in the centroid of 
the line which could also be consistent with a process- 
ing accretion disk. Simulations of emission lines arising 
from an eccentric precessing disk strongly support this 
interpretation (Foulkes et al. 2004). 

In this paper we re-analyze previously published data 
and present new optical spectroscopy of XTE J1118+480 
with the aim of studying and resolving the complex struc- 
ture of the emission line profiles. The observations com- 
prise data acquired at a similar time to those presented 
in Zurita et al. (2002) and a data set acquired 550 days 
after the system entered quiescence. The paper is struc- 
tured as follows: Section 2 presents the observations and 
the data reduction procedure. In Section 3 we determine 
the orbital ephemeris and in Section 4 discuss the spec- 
trum of the donor star. Sections 5 and 6 describe the 
results obtained from the analysis of the Balmer emis- 
sion lines. Finally, in Section 7 we discuss our results 
and a summary is given in Section 8. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed XTE J1118+480 with the Blue Channel 
CCD spectrograph attached to the 6.5 m MMT at the F. 
L. Whipple Observatory. All the spectra were acquired 
with the Loral 3072x1024 CCD binned by two in the 
spatial direction. Other details of the observations are 
summarized in Table 1. 

On six nights during the period 2000 December 1 to 
2001 January 26 UT we used a 500 grooves mm -1 grating 
and a 1.25 arcsec slit width which yielded a dispersion of 
1.18 A pix 1 and a spectral resolution of 3.6 A FWHM. 
On night 2001 March 26 we made use instead of a 300 
grooves mm -1 grating and a 1 arcsec slit which pro- 
vided a dispersion of 1.94 A pix -1 and a resolution of 
5.5 A FWHM. The source visual brightness at the be- 
ginning of the observations was R~18.2 a nd it contin- 
ued t o fade at a rate of 0.003 mag day -1 llZurita et al.1 
2002). Some results obtained from the data acquired 
during the first two nights can be found in McClintock 
et al. (2001b). A final observing run took place on 2003 
January 02/03 at a time when the system was in true 
quiescence (R=18.6). Spectra were acquired using a 832 
grooves mm" 1 grating and a 1.25 arcsec slit. This combi- 
nation gave a dispersion of 0.72 A pix" 1 and a resolution 
of 2.4 A FWHM. 

The images were bias and flat-field corrected with stan- 
dard iraf 3 routines. Spectra were extracted with the 
iraf kpnoslit package. HeNeAr arc lamp exposures 
taken regularly during the observations allowed us to de- 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories. 



rive the pixel-to-wavelength calibration. This was ob- 
tained from cubic spline fits to about 30 arc lines giving 
a root-mean square deviation of typically ^ 0.07 A. The 
wavelength scales of neighbouring arc spectra were inter- 
polated in time and the whole process was checked and 
corrected by using the atmospheric 0[I] 5577.34 A sky 
line, which showed a rms scatter of ^ 9 km s" 1 per night. 
Except on 2003 January 02/03, spectra of flux standards 
were acquired throughout the observations and used to 
correct the target spectra for the instrumental response. 

3. RADIAL VELOCITY CURVE AND ORBITAL EPHEMERIS 

We extracted the radial velocities by cross-correlating 
each target spectrum with archival spectra of six tem- 
plate stars with spectral types ranging through G8V- 
M0V and acquired using the same set-up as for the tar- 
get (see McClintock et al. 2001b). Prior to the cross- 
correlation, the target and template spectra were re- 
sampled into a common logarithmic wavelength scale 
and normalized by dividing with the result of fitting a 
low order spline to the continuum (after masking the 
emission lines in the target spectra). Cross-correlation 
was carried out in the range AA4900 — 6500 and least- 
squares sine fits with all parameters free were performed 
for the radial velocity determinations. Unreliable ve- 
locities obtained from two spectra were not included in 
these fits. We have adopted the K0 V parameters of the 
radial velocity curve because it provided the lowest re- 
duced x 2 (= 3.0) for the fit: 7 = 16 ± 6 km s -1 , K 2 
= 709 ± 7 km s"\ P = 0.1699339 ± 0.0000002 d and 
T = HJD 2451880.1086 ±0.0004, with T corresponding 
to the time of closest approach of the secondary to the 
observer. All quoted uncertainties are 1-er and were ob- 
tained after increasing the error in the radial velocities 
in order to give xt = 1. Figure 1 displays the phase- 
folded radial velocity curve. Our spectroscopic orbital 
period disagrees somewhat (at the 3-c level) with the 
photometric orbital period P ph = 0.169937 ± 0.000001 d 
found by Zurita et al. (2002). Adopting the values of 
Ki and the orbital period, a refined mass function of 
f(M) = 6.3 ± 0.2 M is obtained. Finally, note that by 
using the above refined ephemeris the outburst Doppler 
maps in Torres et al. (2002a) should be rotated clockwise 
by 4° ± 1°, whereas for the tomograms presented in this 
paper the error offsets are negligible for the first two and 
within 2° for the latter. 

4. THE SPECTRUM OF THE DONOR STAR 

Figure 2 presents the spectrum of XTE Jl 118+480 for 
2001 January 26. It was obtained by adding the spec- 
tra with different weights to maximize the signal-to-noise 
ratio after Doppler-correcting them to the rest frame of 
the secondary. By doing this, it is possible to see the 
absorption features from the companion star that would 
otherwise be blurred by the orbital motion of the star. As 
in McClintock et al. (2001b) and Wagner et al. (2001), 
the two strongest photospheric features detected are the 
NaD 1 doublet and the discontinuity in the continuum 
at A5200 caused by the Mgfr triplet, MgH and (proba- 
bly) the TiOa system. Also visible are the MgH band at 
A4780 and the molecular bands of the Ti07 system sit- 
uated longwards of 5840 A. Metallic features such as the 
Ca 1 A4227, 6162, Fe 1 A4384 and the blend of chiefly Fe 
i+Ca 1 at A6495 are also present. The molecular bands 



SPECTROSCOPIC STUDY OF THE BHC XTE J1118+480 



3 



and lines from neutral elements identified are character- 
istic of spectral-types K and M. The spectrum is bluer 
and the absorption features weaker relative to the K5/8 
V template shown in the same figure. This is due to the 
continuum excess from the accretion disk. In this regard, 
the contribution of light from the donor star is ~ 35 ± 8% 
at AA5800-6400 on 2000 December 29 and 2001 January 
26 and 60 ± 10% on 2003 January 02/03. These estima- 
tions were obtained by performing the standard optimal 
subtraction technique (Marsh, Robinson & Wood 1994) 
with K5-M0 dwarf templates. This is comparable to val- 
ues found by other authors. In particular, the relative 
contribution of the secondary was found to be 53 ± 7% on 
2001 April (Zurita et al. 2002) and 45±10% on 2002 Jan- 
uary (McClintock et al. 2003). These two numbers along 
with ours from 2003 January are all consistent with a 
constant ~ 55% contribution from the secondary and in- 
dicate that XTE J1118+480 has reached true quiescence 
on these dates. Finally, there is no clear evidence in the 
night-averaged spectra for absorption from the Lil A6708 
resonance line to an equivalent width (EW) upper limit 
of 0.12 A (1-<t) relative to the observed continuum. In 
comparison, strong Lil has been detected unequivocally 
in the companions of other XRNe. The reader should 
refer to Martin et al. (1994, 1996) for a discussion of the 
possible cause of the anomalously high abundance of Li I 
in XRNe. 

5. THE AVERAGED EMISSION LINE SPECTRUM 

Apart from the late-type companion spectrum, Balmer 
lines up to H<5 in pure emission are detected in our spec- 
tra. Hei and Fen emission lines which are commonly 
observed in XRNe and DNe in quiescence are not dis- 
cernible, although it has been shown that at least Hei 
A5876 is present in the spectrum of XTE Jl 118+480 but 
veiled by the continuum from the companion star (Mc- 
Clintock et al. 2003). In order to study the Balmer 
line profiles, we averaged the spectra in the observer's 
rest-frame for the nights when 6 or more spectra were 
acquired. The line profiles have a FWHM of ~ 2500 km 
s _1 and appear double-peaked, a profile shape which 
is indicative of emission arising from an accretion disk. 
Table 2 presents the measured full width zero intensity 
(FWZI) and EWs for the lines. In Table 3 we give the 
position of the peaks for Ha obtained by fitting the av- 
eraged line profiles with a 3-Gaussian function, with one 
used to account for the base of the line and the other 
two for the narrow peaks in the line. Also their dif- 
ference (peak-to-peak separation) and mean (centroid of 
the line) are listed. From Table 2 it is obvious that 
the EWs measured for Ha during 2003 January 02/03 
are higher than during approach to quiescence. This is 
caused in part by the lower contribution from the ac- 
cretion disk to the optical continuum during quiescence 
(Section 4). By looking at Table 3, there is no evidence 
for a change in the peak-to-peak separation during the 
period of our observations, having a weighted mean value 
of 1780 km s _1 and a rms scatter of 60 km s _1 . In- 
stead, during the approach of XTE J1118+480 to qui- 
escence, the position of the centroid of the line with re- 
spect to the rest wavelength varies noticeably from night 
to night with an amplitude of 250 km s _1 . The centroid 
also shows displacements during quiescence. Addition- 
ally, the averaged Ha profiles during the decline to qui- 



escence show long-term changes in the relative intensity 
of the red and blue peaks, a fact previously reported by 
Zurita et al. (2002) and also observed during outburst 
(Torres et al. 2002a). During quiescence the peaks of the 
line are slightly asymmetric. 

6. THE TIME-RESOLVED Ho LINE PROFILE 

In order to study the orbital variations in the Ha emis- 
sion line profile, we firstly determined the EWs for the in- 
dividual spectra on those nights with good orbital phase 
coverage. The resultant Ha EWs are plotted as a func- 
tion of the orbital phase in Figure 1. During 2001 Jan- 
uary 26, the EW reaches maximum at orbital phase 0.1 
and no obvious minimum is observed. During 2003 Jan- 
uary 02/03, the Ha EWs display a maximum at phase 
0.5. The plot also shows that the EW is significantly 
higher around phase 0.5 for the first night in compari- 
son to the second night. The variations are unlikely to 
be caused by the ellipsoidal modulation in the contin- 
uum of the secondary star since one would expect to see 
two equal minima in the EWs at phases 0.25 and 0.75 
corresponding to the two maxima in the ellipsoidal mod- 
ulation. The variations more probably represent intrinsic 
changes in the Ha line flux with orbital phase. 

The upper panels in Figure 3 show the evolution of 
the Ha profile with orbital phase in the form of trailed 
spectrograms. They also illustrate the prominent long- 
term changes in the velocity shifts relative to the rest 
wavelength of the line and in the asymmetry of the peak 
intensities (Section 5). Looking to the panels, it is pos- 
sible to trace the narrow emission components present 
in the double-peaked line profiles. During 2000 Decem- 
ber 29, a sharp emission component is visible at orbital 
phases 0.3-0.7 which crosses from red to blue at phase 
~ 0.4. On 2001 January 26 the line profile shows an 
emission component in the phase interval 0.7-0.9 which 
crosses from blue to red at phase ~ 0.9. During 2003 
January 02/03 an emission component crosses zero ve- 
locity from red to blue at around phase 0.5. Because 
of the gap around phase 0.0 in the orbital phase cover- 
age for the latter epoch, it is not possible to know if the 
narrow component is observable at all orbital phases. 

Since the narrow components in the emission line pro- 
files are weak and show a complex behaviour, a multi- 
Gaussian fit to each individual spectrum will render poor 
information about their origin. However, we can gain in- 
sight into the Ha emitting regions in XTE J1118+480 
by using the Doppler tomography technique (Marsh & 
Home 1988). This technique reconstructs the brightness 
distribution of the binary system in velocity space, al- 
lowing us to localize the emission structures which are 
not easily recognizable in the individual spectra. We use 
the maximum-entropy method (MEM) of building the 
tomograms on those nights with orbital phase covering 
at least half an orbital cycle and with good spectral res- 
olution: 2000 December 29, 2001 January 26 and 2003 
January 02/03. The bottom panels in Figure 3 display 
the derived MEM tomograms with the theoretical path 
of the gas stream and the Keplerian velocities of the disk 
along the stream for a radial velocity semi-amplitude of 
the companion star of K 2 =709 km s _1 (Section 3) and a 
mass ratio q = 0.037 (Orosz 2001). 

During the approach to quiescence, enhanced emission 
in the -Vx, +Vy (upper-left) quadrant together with 
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the emission arising from the rotating accretion disk is 
easily recognizable in the Doppler images. The loca- 
tion of the former emission, halfway between the the- 
oretical disk and gas stream kinematic paths, is sim- 
ilar to that observed in the Balmcr and Hen A4686 
maps during outburst (Torres et al. 2002a). However, 
now there is a greater overlap between this region and 
the theoretical trajectories. On 2000 December 29, a 
bright spot is present between both theoretical paths at 
0.5 ± 0.1 Rlj = 0.4 ± 0.1 a where Rlj is the distance of 
the compact object from the inner Lagrangian point and 
a the binary separation. The position was determined 
by using lines passing through points in the theoretical 
paths at equal distance from L\. On 2001 January 26 a 
bright spot is again present, but now it is approximately 
placed on the ballistic trajectory of the gas stream at 
0.77 ± 0.06 R Ll = 0.61 ± 0.05 a. Finally, note the en- 
hanced emission in the + Vx , -Vy (right-lower) quadrant 
for both tomograms. Its origin will be discussed in the 
next section. 

The tomogram computed during quiescence shows no 
clear evidence of emission from a bright spot on the ring- 
like distribution of the accretion disk. An emission com- 
ponent crossing the double-peaked profile is evident in 
the trailed spectra and therefore it must be reflected in 
the tomogram at a position defined by the amplitude and 
phase of its modulation. In order to find its counterpart 
in the tomogram, we subtracted the azimuthally sym- 
metric part of the map. This was obtained by taking 
the median of the image values within annuli centered 
on the center of symmetry of the ring emission, which 
falls at - (-100,0) km s" 1 rather than (0,-#i). Fig- 
ure 4 shows the result of the subtraction. We can see 
that the core of the residual emission is situated near the 
predicted Roche lobe of the secondary. The spread is 
approximately consistent with the combined effect of the 
finite spectral resolution (110 km s _1 ) and orbital phase 
smearing (ranging 20 to 270 km s _1 ). 

7. DISCUSSION 

From spectroscopic observations of XTE Jl 118+480 
during its approach to quiescence, Zurita et al. (2002) 
found that the variability in the centroid of the nightly 
averaged Ha profile was consistent with a periodic mod- 
ulation of 26 days (half of the 52 day precession period 
for the accretion disk implied by the superhump period) , 
although a 52 day periodicity was not ruled out. Combin- 
ing our centroid measurements obtained for Ha during 
approach to quiescence (Table 3) with those of Zurita et 
al., a modulation on the precession period is definitely 
favoured and the 26 day periodicity is ruled out (Fig- 
ure 5). The centroids measured for 2003 January 02/03 
(Table 3) are also displaced with respect to the systemic 
velocity which suggests that the accretion disk contin- 
ues precessing during true quiescence. This is possible 
only if the outer disk radius is large enough to reach or 
exceed the 3:2 resonance tidal instability. By using the 
Ha mean pcak-to-peak separation (Aw) and the semi- 
amplitude of the secondary, we estimate an outer radius 
of R d /a = (1 + q) {2K 2 /Avf = 0.66 ± 0.05 (see e.g. 
Dubus et al 2001). This is a value definitely larger than 
the 3:2 resonance radius 132/ a (l) = 0.47 (Whitchurst & 
King 1991) and comparable to the expected disk tidal 
truncation radius Rx/a = 0.6/(1 + q) = 0.58 (Warner 



1995). Additionally, the hotspot in the 2001 January 26 
tomogram is located at a position consistent with the 
predicted outer disk radius. 

Before continuing with the discussion on the Doppler 
maps obtained in the previous section, it is important to 
stress that the Doppler tomography technique assumes 
(in particular) that all the emitting regions co-rotate in 
the binary and that the flux from any of these regions is 
constant with time. These assumptions might not be ful- 
filled in a precessing eccentric disk during the course of a 
single orbit (Foulkes et al. 2004) and therefore the results 
obtained from the tomograms must be interpreted with 
caution. Bearing this in mind, the tomograms for 2000 
December 29 and 2001 January 26 show emission from 
the hotspot and the gas stream illustrating that mass 
transfer from the donor star has continued during the 
approach of the system to quiescence. They also show 
changes in the position of the hotspot (see Section 6). 
It is interesting to note here that the changes observed 
in the position of the hotspot may be modulated with 
the precession period because of the varying geometry of 
the elliptic disk. On the other hand, the emission in the 
lower-right quadrant may be interpreted as arising from 
an overflowing gas stream which after the first impact in- 
tersects again the disk close to phase 0.5 (see e.g. Kunze, 
Speith & Hessman 2001) or as arising from a non- uniform 
brightness distribution across the disk. In the latter case, 
the formation of a region of enhanced density (tempera- 
ture) in the disk may be caused by tidal forces from the 
companion star (see e.g. Steeghs & Stehle 1999). It is 
clear from the long-term variations in the centroid and 
peak intensities of the Ha profile that the disk in XTE 
Jl 118+480 continues to experience the tidal stresses from 
the companion star. 

The residual Ha image obtained by subtracting the 
symmetric part from the 2003 January 02/03 Doppler 
image uncovers emission concentrated towards the ex- 
pected position of the secondary star. Ha emission from 
the hotspot or gas stream is not detected, a fact that 
points to a decrease or pause in the mass transfered from 
the donor star to the accretion disk. A similar scenario 
was suggested to explain the quiescent Ha tomograms of 
Nova Muscae 1991 (Casares et al. 1997) and Centaurus 
X-4 (Torres et al. 2002b) which show the same pattern 
(secondary's emission present/ absence of hotspot). In 
order to constrain the nature of the Ha emission from 
the secondary star, we isolated the narrow line compo- 
nent of the line profile as follows: first, an accretion disc 
profile was created by averaging all the target spectra in 
the observer's rest frame. Next, seventeen copies of the 
simulated disc profile were shifted to the rest frame of the 
secondary star at the times of the individual target spec- 
tra and then averaged using the same weights as in the 
first step of the procedure. Finally, the resultant profile 
was subtracted from the Dopplcr-corrccted Ha averaged 
profile of XTE Jl 118+480. This left us with a residual 
with an EW=1.0 ± 0.5 A which represents 1% of the 
total Ha emission. For a R quiescent magnitude of 18.6, 
a distance of 1.8 kpc and assuming a secondary with a 
stellar radius equal to the volume radius of its Roche lobe 
(R 2 = R h (2) = 0.15 a, with a = 2.5 R for Mi = 7 M Q 
and i — 80°), we derive from the EW of the Ha residual 
a Ha flux at the surface of the star (Fh q ) of 3 x 10 6 erg 

— 9 — 1 

cm z s . 
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Docs X-ray irradiation power the Ha emission from 
the donor star? The intrinsic quiescent X-ray luminos- 
ity in the 0.3-7 keV band was estimated by McClintock 
et al. (2003) to be L x « 3.5 x 10 30 erg s^ 1 . The 
implied X-ray irradiation at the secondary star is then 
Fx = Lx/i^a 2 ) = 9x 10 6 erg cm~ 2 s -1 which can ac- 
count for the Ha emission if at least a 30% of the incident 
X-ray flux is reprocessed to Ha photons. We can obtain 
a rough estimate of the likely fraction of the X-ray lu- 
minosity converted to Ha photons by assuming that the 
Ha emission line observed in X-ray bright LMXBs is due 
completely to X-ray reprocessing in the disk. Then: 



EW(A) 



(7000 - 3000 A) 



x L. 



opt 



10- 



J X 



where we have taken into account that for LMXBs 
the averaged optical to X-ray luminosity is L opt (3000 — 
7000 A)/L x (2 - 11 keV) - 0.002 (van Paradijs & Mc- 
Clintock 1995) and the EW of the Ha emission line is of 
the order of a few A (see e.g. Shahbaz et al. 1996; Torres 
et al. 2002a). Therefore we find it improbable that the 
X-ray irradiation could be the source of the Ha emission 
in the secondary star. 

A chromospherically active secondary is an alterna- 
tive explanation to the Ha emission. We can com- 
pare the above estimations with the results reported by 
Soderblom et al. (1993) for rapidly rotating dwarfs in the 
Pleiades. In order to compare the EW of the Ha resid- 
ual, we have first subtracted the underlying Ha profile 
of the secondary star from the Ha residual by using the 
dwarf templates and taking into account the veiling from 
the accretion disc (Section 4). This operation yielded a 
corrected EW of 3 A. By looking at table 6 of Soderblom 
et al. we find that the Ha emission in the companion of 
XTE J1118+480 is comparable in strength to that ob- 
served in some of the K dwarfs in their sample. For ex- 
ample, for the K7V star listed as k403, EW=3.41 A and 



6.12 x 10° erg cm 1 s 



7.1. Comparison of XTE J1118+480 and GRO 
J0422+32 

Among the black hole XRNe, GRO J0422+32 (Nova 
Per 1992; P=5.1 hr) has the closest orbital period to 
XTE J1118+480 and has also exhibited superhumps dur- 
ing the 1992 outburst (O'Donoghue & Charles 1996 and 
references therein), suggesting a precessing accretion disk 
in the system in accordance with its low mass ratio 
(q = 0.12±°;°^; Harlaftis et al. 1999). In contrast, the de- 
cay to quiescence was interrupted by two mini-outbursts 
in 1993 August and December. Motivated by the similar- 
ities between these systems, we have searched the litera- 
ture for evidence of a precessing accretion disk from the 
Ha emission line of GRO J0422+32. We have found that 
asymmetries in the double-peaks of the averaged line pro- 
file have been observed (see Table 4) which we suggest 
are most likely due to a precessing accretion disk. Addi- 
tionally, the centroid of the emission profiles should also 
vary if the disk is precessing. This may be reflected in 
the Ha systemic velocity of +142 ±4 km s _1 obtained by 
Garcia et al. (1996), in comparison with the secondary's 
systemic velocity of 7 = 11 ± 8 km s _1 (Harlaftis et 
al. 1999). Note that Filippenko, Matheson & Ho (1995) 



found a Ha systemic velocity of 26.0 ± 1.6 km s _1 by 
combining data obtained on 1994 November 8/9, 1995 
January 26/27. Qualitatively, the Ha and Hen A4686 
trailed profiles in fig. 6 of Casares et al. (1995) are 
red-shifted by ~ 100 km s _1 . We suggest that these 
variations are also associated with a precessing disk. Fu- 
ture spectroscopic monitoring of XRNe in outburst and 
decline will enable us to confirm unambiguously the rela- 
tion between the long-term variations in the line profiles 
and the precession of the accretion disk. 

8. SUMMARY 

We have presented comprehensive optical spectroscopy 
of the XRN XTE Jl 118+480 from 2000 December to 
2003 January. From the analysis of the absorption 
line spectrum, we have established an orbital period of 
P = 0.1699339 ± 0.0000002 d and a radial velocity semi- 
amplitude of the secondary of if 2 = 709 ± 7 km s _1 . The 
implied updated mass function is /(M) = 6.3 ± 0.2 M Q . 
Inspection of the averaged spectrum of XTE J1118+480 
in the rest frame of the secondary star exhibits numerous 
photosphcric features typical of a secondary of spectral 
type K-M. Our analysis of the absorption features sug- 
gests that the system has reached true quiescence with 
the companion star contributing <~ 55% of the total flux 
over the range AA5800 — 6400. The long-term variations 
in the centroid of the Ha profile during the decline to 
quiescence show a periodicity of <~ 52 days and an am- 
plitude of 250 km s _1 . They are interpreted as due to a 
precessing accretion disk. Comparison between two to- 
mograms of the Ha line profile during the approach of 
the system to quiescence show a hotspot whose position 
varies. The change in the hotspot location may be re- 
lated to the changes in the gas stream impact site at the 
edge of the disk as it precesses. During quiescence, both 
hotspot and gas stream are not present which points to 
a decrease or cessation of the mass transfer from the sec- 
ondary. Ha emission is detected near the location of the 
secondary which is most probably due to chromospheric 
activity. 

We thank Jeff McClintock who generously made avail- 
able his XTE J1118+480 data. This work was sup- 
ported in part by NASA Contract NAS8-39073 and LISA 
Grant NAG5- 10889. Use of MOLLY, DOPPLER and 
TRAILER routines developed largely by T. R. Marsh is 
acknowledged. We are grateful to the anonymous referee 
for useful comments which improved the quality of the 
manuscript. 
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Fig. 1. — Upper panel: a total of 85 radial velocities of the secondary star in XTE J1118+480 folded on the ephemcris of Sect. 3.1. 
The best sine- wave fit is also shown. Lower panel: EWs of Ha versus orbital phase for 2001 January 26 (solid triangles), 2003 January 02 
(squares) and 2003 January 03 (solid squares) 
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Fig. 2.— Average optical MMT spectrum of XTE J1118+480 on 2001 January 26 (UT) together with the spectrum of a K8/5 V template. 
The spectra have been normalized at A4100 and the spectrum of XTE J1118+480 has been shifted vertically by +2 for the sake of clarity. 
Major features are identified. The Earth symbol denotes the atmospheric features and bl the metallic blend. 

2000 Dec 29 2001 Jan 26 2003 Jan 02,03 




-1000 1000 -1000 1000 -1000 1000 



Velocity (km/s) Velocity (km/s) Velocity (km/s) 

Fig. 3. — Top panels present the trailed spectrograms of Ha. Empty strips represent gaps in the phase coverage. For the sake of clarity, 
for 2000 December 29 and 2001 January 26 the same cycle has been plotted twice. For 2003 January, the lower and upper half of the 
spectrogram arc data from January 02 and 03 respectively. Lower panels show the computed MEM Doppler maps. For 2003 January 02/03, 
we show the tomogram obtained by combining the data acquired during both nights. The Roche lobe of the secondary star, the predicted 
velocities of the gas stream (lower curve) and of a Keplerian disk along the stream (upper curve) are plotted. Distances in multiples of 
O.lRii arc marked along both curves with open circles. The center of mass of the system is denoted by a cross. 



SPECTROSCOPIC STUDY OF THE BHC XTE J1118+480 



9 




Fig. 5. — The changing centroid of the Hct emission line during approach to quiescence. The velocity points without errors have been 
extracted from fig. 5 of Zurita et al. The best sine-wave fit (P~ 52 d) is also shown 
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TABLE 1 

Journal of Observations 



Date 


No. 


Exp. time 


HJD start 


HJD end 


A range 


(UT) 


spectra 


(v 


/ I O A ET 1 c\r\c\ \ 

(+2,451,000.) 


(+2,451,000.) 


( A ) 


2000 Dec 1 


2 


900 


880.0097 


880.0224 


3650-6825 


2000 „ 4 


6 


900-1200 


882.5471 


882.9630 


3940-7120 


2000 „ 24 


6 


1200 


902.8517 


902.9851 


4100-7015 


2000 „ 28 


3 


1100-1200 


907.0248 


907.0543 


4090-7005 


2000 „ 29 


13 


1200 


907.8632 


908.0551 


4085-7000 


2001 Jan 26 


38 


600 


935.7534 


936.0420 


4075-6990 


2001 Mar 26 


2 


1200 


994.7991 


994.8147 


3715-8510 


2003 Jan 02 


8 


900 


1,641.9798 


1,642.0635 


4925-6850 


2003 Jan 03 


9 


900 


1,642.9625 


1,643.0581 


4930-6860 



TABLE 2 

FWZI a and Equivalent Widths 3 - for the Balmer lines 





Ho 


H/3 


H 7 




UT Date 


FWZI 


EW 


FWZI 


EW 


FWZI 


EW 


FWZI EW 


2000 Dec 4 


88±5 


-22+2 


46+10 


-8+2 


78+5 


-13+5 


70+5 -12+10 


2000 Dec 24 


78±5 


-27+5 


55+6 


-12+3 


60+5 


-9+3 




2000 Dec 29 


85+6 


-26+4 


65+10 


-7+3 


76+5 


-16+5 




2001 Jan 26 


85+5 


-32+6 


51+4 


-11+4 


76+5 


-16+9 




2003 Jan 02 


90+2 


-76+8 












2003 Jan 03 


99+2 


-70+6 













Note. — The EWs are the mean of the values obtained for the individual profiles and the uncertainties correspond to the standard 
deviation. The uncertainties in the FWZIs were estimated by looking at the scatter in the values when selecting different wavelength 
intervals to set the local continuum level. 

a In angstroms. 



TABLE 3 
Fits to the Ha Line Profile 



Date 


v 6 


Vr 


V r -V 6 


(V b +V r )/2 


UT 


(km s ) 


(km s 1 ) (km s ) 


(km s ) 


2000 Dec 4 


-1163+34 


664+14 


1827+37 


-250 ±18 


2000 Dec 24 


-644+16 


1021+23 


1665+28 


188 ±14 


2000 Dec 29 


-651 ±8 


1120+17 


1771+19 


234 ±9 


2001 Jan 26 


-1134+12 


627+7 


1761+14 


-254 ±7 


2003 Jan 02 


-952+6 


811+8 


1763+10 


-71 ±5 


2003 Jan 03 


-1009+6 


823+8 


1832+10 


-93 ±5 



Note. — Vt and V r designate the positions (shifts) respect to the rest wavelength of the line of the blue and red peaks respectively. 
The uncertainties were calculated after scaling the error bars to give xt = 1- 
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TABLE 4 

Asymmetries in the averaged double-peaked Ha profile of GRO J0422+32 



Night/s 


Net Exposure Stronger Comments 


refs. 


averaged 


(hr) 


peak 






1993 Oct 10 


5.3 


Blue 


Decay to quiescence after August mini-outburst 


1 


1993 Oct 12 


5.3 


Blue 






1993 Dec 9 


6 


Red 


Maximum December mini-outburst 


2 


1993 Dec 10 


5 


Red 






1993 Dec 16-19 


~ 11 


Blue 


Initial decline of mini-outburst 


3 


1994 Nov 5-7 


> 14.5 


Blue 


Quiescence 


4 


1995 Jan 26 


4.9 


Red 


Quiescence 


5 



Note. — References: 1.- Garcia ct al. (1996); 2.- fig. 5d in Callanan et al. (1995): 3.- fig. 2 of Casarcs ct al. (1995); 4.- fig. 2 in Orosz 
& Bailyn (1995): 5.- fig. 1 in Harlaftis ct al. (1999) 
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